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ABSTRACT 


The  M3 3  projectile  has  been  analyzed  using  the  ANSYS  CFX  code  that  is  based  on  the 
numerical  solution  of  the  full  Navier-Stokes  equations.  Simulation  data  were  obtained 
against  various  Mach  numbers  ranging  from  M=  0.5  to  M=  2.6  at  0°  and  2°  angles  of 
attack.  Simulation  data  were  also  obtained  against  various  angles  of  attack  from  0°  to  85° 
for  M=0.5. 

For  Mach  numbers  between  M=  0.5  to  2.6,  the  results  obtained  using  the 
combined  k-epsilon  and  Shear  Stress  Transport  model  show  good  agreement  with  the 
experimental  range  data  for  the  normal  force  and  pitching  moment  coefficient.  The  drag 
coefficient  at  zero  angle  of  attack  tended  to  be  over  predicted  by  an  average  error  of 
1 1.6%  with  the  highest  error  occurring  at  M=  1.5. 

For  varying  angle  of  attack  up  to  85°  at  M=  0.5,  the  results  obtained  from  CFX 
code  were  compared  with  simulation  results  obtained  from  AP09.  The  data  showed  good 
agreement  only  up  to  20°  angles  of  attack. 
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I. 


INTRODUCTION 


The  employment  of  artillery  to  strike  strategic  targets  or  simply  to  provide 
overwhelming  fire  power  over  adversaries  has  been  one  of  the  key  tactics  employed  in 
the  history  of  combat.  In  view  of  this,  the  design  of  projectiles  achieving  longer  ranges 
and  better  accuracies  was  viewed  as  the  main  thrust  of  this  research  and  development.  To 
achieve  this,  we  need  to  understand  the  effects  on  the  aerodynamic  properties  of  the 
projectile  during  its  flight  with  specific  interest  against  the  variation  of  speed  and  angle  of 
attack. 

Conventional  approach  to  predict  the  aerodynamic  properties  of  a  projectile  is 
through  wind  tunnel  and  actual  range  testing.  Alternate  methods  using  Computational 
Fluid  Dynamic  (CFD)  software  or  other  computational  codes  are  also  able  to  predict 
these  properties  fairly  well.  This  report  aims  to  validate  the  aerodynamic  properties  of  the 
M3 3  projectile  at  different  Mach  numbers  as  well  as  at  various  angles  of  attack  using  the 
CFX  code. 

The  ANSYS-CFX  code  is  a  commercial  CFD  program  used  to  simulate  fluid  flow 
in  a  variety  of  applications  such  as  gas  turbine  engines  and  aircraft  aerodynamics.  The 
numerical  solution  techniques  employed  are  based  on  the  Reynolds-averaged  Navier- 
Stokes  equations  and  the  Euler  equations.  Results  obtained  through  the  CFX  code  will 
also  be  compared  with  previous  results  obtained  by  other  codes  such  as  the 
Aeroprediction  Code  2009  (AP09)  and  the  Parabolized  Navier-Stokes  (PNS)  viscous 
flow  solver.  The  results  will  also  be  compared  with  experimental  range  data  to  assess  the 
codes  accuracy. 

In  view  of  the  extensiveness  of  the  research  requirements,  this  report  will  only 
study  the  aerodynamic  properties  of  a  static  M3  3  in  flight.  The  effects  of  a  rotation  in  the 
case  of  a  spin-stabilized  projectile  will  not  be  studied.  Due  to  the  unclassified  nature  of 
the  study,  research  will  also  be  limited  to  only  unclassified  information  or  information 
from  open-source  documents. 
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II.  LITERATURE  REVIEW 


Conventional  approaches  of  aerodynamic  prediction  on  flight  bodies  are  done 
using  wind  tunnel  and  actual  range  testing.  This  approach  is  both  time  consuming  and 
cost  intensive.  With  the  development  of  modem  CFD  methods  and  software,  engineers 
are  able  to  achieve  significant  accuracies  in  the  prediction  of  aerodynamic  properties 
without  the  hassle  of  time  and  cost. 

There  are  numerous  unclassified  studies  and  researches  done  on  the  aerodynamic 
prediction  of  projectiles.  The  Army  Research  Laboratory  (ARL)  provides  a  wide  range  of 
experimental  data  and  simulation  results  on  such  studies.  Aerodynamic  predictions  using 
traditional  codes  such  as  MCDRAG  and  NSWCAP  for  drag  predictions  of  projectiles  at 
transonic  and  supersonic  speeds  have  already  been  done  since  the  1980s  [1].  More 
advanced  computational  methods  such  as  PNS  and  AP09  have  also  achieved  fairly 
accurate  results  [2]  -  [3].  This  provides  a  strong  baseline  for  comparison  of  the  results 
obtained  with  the  CFX  code. 

With  the  development  of  the  ANYSYS-CFX  code,  engineers  can  also  study  the 
flow  field  and  its  relation  to  the  aerodynamic  properties  of  the  projectile.  It  is  the 
objective  of  this  report  to  provide  data  on  the  aerodynamic  properties  of  the  M33 
projectile  at  various  speeds  and  angles  of  attack  using  the  CFX  code.  Solution  techniques 
and  procedures  for  this  report  are  referenced  to  the  previous  study  done  by  Mr  Tan  Wei 
Chieh  [4].  Results  from  this  report  will  be  matched  against  experimental  data  available 
for  the  M3 3  projectile  found  in  [5],  [6]  and  [7],  as  well  as  the  computational  data 
obtained  from  the  AP09  program.  The  processes  and  procedures  in  the  development  of 
this  report  will  also  be  documented  in  detail  for  future  course  works  and  analytical 
studies  using  the  CFX  code. 
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III.  SOLUTION  TECHNIQUES 


A.  THEORETICAL  APPROACH 

1.  General 

For  flow  over  bodies  of  revolution,  the  slender  body  theory  is  applicable  to  yield 
reasonable  estimates  of  the  aerodynamic  forces  and  moments  provided  that  the  small 
perturbation  theory  is  satisfied  [8]. 

2.  Governing  Equations 

F,  =  pJ^7ra[R\^,^  -  R\^J  (1) 

Using  the  base  area  as  the  reference  area,  Equation  (1)  leads  to 

C,=la  (2) 

It  is  understood  that  no  lift  will  be  generated  at  zero  angle  of  attack.  The  equations 
also  imply  that  for  a  body  of  constant  diameter  or  a  body  with  a  pointed  nose  and  tail,  no 
lift  will  be  generated  as  well  [8].  The  slender  body  theory  provides  rough  estimates  for 
the  aerodynamic  properties  but  more  accurate  analysis  are  usually  required  in  practice. 

3.  Conical  Shock  Wave 

Since  lift  is  primarily  generated  by  the  nose  cone,  it  is  imperative  to  understand 
the  aerodynamic  properties  occurring  at  the  nose  during  flight.  The  study  of  the  conical 
shock  wave  has  been  done  by  Jones  [9]  and  Kopal  [10]  to  satisfy  the  non-linear 
differential  equations  between  the  shock  and  the  body  at  zero  incidences  to  the  free- 
stream.  The  solutions  were  later  computed  by  Sims  [1 1]  in  1964. 

B.  COMPUTATIONAL  APPROACH 

1.  General 

The  CFX  simulation  was  done  using  a  half-bodied  symmetry  model.  This  was 
achieved  by  cutting  out  a  half-bodied  M33  projectile  configuration  from  the  control 
volume.  This  approach  will  show  the  flow  field  across  the  projectile  as  well  as  save 
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computational  time.  Flow  is  simulated  by  setting  the  front  of  the  control  volume  facing 
the  tip  of  the  projectile  as  an  inlet. 


2.  Governing  Equations 


The  governing  equations  for  the  CFX  code  are  the  conservation  equations  for 
mass,  momentum  and  energy  for  a  compressible  medium.  Equation  (3)  is  the  mass 
conservation  equation  (continuity  equation)  and  Equations  (4)  are  the  momentum 
conservation  equations  (Navier-Stokes  equations)  for  compressible  flow.  Together  with 
the  energy  equation  in  Equation  (5)  and  the  equation  of  state  in  Equation  (6),  these  six 
equations  constitute  a  system  to  solve  for  the  six  unknowns  that  are  the  velocity 
component  u,  v,  w,  the  pressure,  the  density  and  the  temperature. 
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3.  Boundary  Conditions 

The  surface  of  the  projectile  is  set  as  an  adiabatic  wall  under  no-slip  conditions. 
The  inlet  and  outlet  boundary  condition  are  set  as  free-stream  conditions  while  a 
symmetry  condition  is  applied  at  the  centerline  axis  of  the  projectile.  The  rest  of  the 
control  surfaces  are  set  as  opening.  The  default  flow  field  is  set  as  ffee-stream  conditions. 

4.  Turbulence  Model 

The  k-epsilon  (k-  s)  and  Shear  Stress  Transport  (SST)  model  were  explored  for 
the  turbulence  modeling  of  the  simulation. 

The  k-  s  model  is  a  two  equation  model  using  the  turbulent  kinetic  energy,  k  and 
the  turbulent  dissipation,  s.  It  accounts  for  history  effects  like  convection  and  diffusion  of 
turbulent  energy  that  makes  it  a  reasonably  accurate  prediction  model  although  it  does 
not  perform  well  in  cases  with  large  adverse  pressure  gradients  and  flow  separation  [12]. 

The  SST  model  combines  the  effects  of  the  k-omega  (k-co)  model  which  uses 

specific  dissipation,  co  and  the  k-  s  model.  The  k-co  model  is  used  in  the  inner  parts  of  the 

boundary  layer  while  the  k-  s  model  is  used  in  the  free-stream  as  the  prior  model  is  too 

sensitive  to  the  inlet  ffee-stream  turbulence  properties.  The  SST  model  performs  well  for 
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cases  of  adverse  pressure  gradients  and  separated  flow  but  tends  to  over-predict 
turbulence  level  in  regions  with  large  normal  strain  and  strong  acceleration  [12], 

5.  Free-stream  Flow 

The  free-stream  velocity  is  simulated  as  the  u-component  of  the  inlet  boundary 
condition  for  zero  angle  of  attack.  At  other  angle  of  attack  with  reference  to  the  X-Y 
plane,  the  inlet  flow  is  divided  into  the  u-component  and  v-component  as  shown  in 
Figure  1.  This  approach  will  simulate  flow  at  different  angles  of  attack  and  save 
computational  time. 


The  u  and  v  component  were  derive  from  the  free-stream  Mach  number  using 
Equation  (7). 

u  =  Mach  X 340.3  x  cos  a 

V  =  Mach  X 340.3  x  sin  a  (7) 

6.  Results 

Fy  and  Fx  obtained  from  the  CFX  simulations  represents  the  axial  and  normal 
force  respectively.  At  zero  angle  of  attack,  Fy  and  Fx  are  equivalent  to  the  lift  (Fl)  and 
drag  force  (Fd)  respectively.  For  other  angles  of  attack,  Fy  and  Fx  are  not  equal  to  Fl  and 
Fd  as  shown  in  Figure  2.  Fl  and  Fd  will  have  to  be  obtained  through  further  calculations 
as  shown  in  Equation  (8). 
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Figure  2.  Force  Diagram  (CFX). 


Fjy  =  cos  a+  Fy  sin  a 
F^  =  Fy  cos  cc-  F^  sin  a 


(8) 


After  obtaining  Fy,  Fx  and  the  torque  with  respect  to  z-axis  (Tz),  the  results  for  Cl 
and  Cd  at  various  angles  of  attack  were  derived  using  Equation  (9)  while  Cno,  Cjvia  at 
various  Mach  numbers  were  derived  using  Equation  (10),  (1 1)  and  (12).  Cdo  is  simply  the 
drag  coefficient  at  zero  angle  of  attack.  0°  and  2°  were  used  for  the  computation  for  CNa 
and  Cmu-  Note  that  all  the  forces  obtained  from  the  simulation  need  to  be  multiplied  by 
two  as  the  model  only  accounts  for  a  half-bodied  projectile.  The  length  from  the 
projectile  base  to  the  center  of  gravity  (0.023m)  was  used  as  the  reference  length  for  the 
computation  of  Cm- 


(9) 


cos  a  +  C^  sin  a 


(10) 
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IV.  PROCESS  AND  PROCEDURES 


A.  SOLIDWORKS  2010 

1.  Introduction 

SolidWorks  2010  is  a  3D  mechanical  CAD  software  used  for  design  and 
modeling  purposes.  It  was  released  in  December  2010  and  has  been  upgraded  to  its  latest 
version  of  SolidWorks  2012  in  September  2011.  It  uses  a  parametric  featured-based 
method  to  create  models  and  assemblies.  Details  of  the  processes  and  input  parameters 
can  be  found  in  Appendix  A. 

2.  Two  Dimensional  Model 

The  2D  model  of  the  M33  projectile  as  shown  in  Figure  3  was  drawn  with 
reference  to  the  geometry  inputs  in  [6]. 


The  ogive  nose  profile  of  the  M3 3  projectile  was  drawn  using  tangent  arcs  for  the 
nose  cone  and  nose  tip.  The  control  volume  was  set  as  10  times  the  reference  length  of 
the  projectile  to  avoid  wall  interference  effects  as  is  often  seen  in  wind  tunnel  testing. 
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3. 


Three  Dimensional  Model 


There  are  two  ways  to  produce  the  3D  model  of  the  M33  projectile  for  simulation. 
One  way  is  to  generate  the  projectile  as  a  full  solid  3D  model  as  shown  in  Figure  4  and 
the  other  is  to  cut  out  the  half-body  of  the  projectile  from  the  control  volume  as  a  cavity 
as  shown  in  Figure  5.  The  latter  method  was  chosen  for  simplicity  and  saving  of 
computational  time. 


Figure  4.  Solid  3D  Model  of  M33  Projectile  (SolidWorks). 


The  model  was  saved  as  a  parasolid  file  to  be  imported  into  the  geometry 
interface  of  the  ANSYS-CFX  code. 

B.  CFX 

1.  Introduction 

The  CFX  code  uses  the  full  Navier-Stokes  equations  to  solve  for  fluid  flows  over 
a  region  of  interest  with  given  specific  boundary  conditions.  It  complements  traditional 
empirical  methods  and  wind  tunnel  testing  as  an  alternate,  cost-effective  tool  with  vivid 
graphic  interface  for  the  study  of  fluid  flow.  Details  for  the  processes  and  input 
parameters  can  be  found  in  Appendix  B. 

2.  Mesh 


A  portion  of  the  mesh  for  over  650,000  nodes  and  2  million  elements  is  seen  in 
Figure  6. 
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Figure  6.  Mesh  of  M33  Projectile  (CFX). 


The  parameters  for  the  mesh  input  differ  from  model  to  model.  Engineers  need  to 
discover  the  optimal  input  parameters  to  achieve  high  fidelity  through  trials  and  errors. 
The  fundamentals  for  better  fidelity  were  to  reduce  the  element  size  and  achieving  a 
small  Y-plus  value  after  post-processing.  Reducing  the  element  size  refines  the  mesh 
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especially  in  areas  with  protruding  angles.  This  will  reduce  the  computational  error.  Y- 
plus  is  the  non-dimensional  wall  distance  for  a  wall-bounded  flow  that  is  commonly  used 
in  boundary  layer  theory  and  in  defining  the  law  of  the  wall  [13]. 

To  reduce  element  size,  engineers  can  try  to  reduce  the  sizing  input  parameters  for 
minimum  and  maximum  size  along  with  the  curvature  angle  in  the  mesh  module.  This 
will  produce  a  model  with  higher  number  of  nodes  and  elements  but  need  not  necessarily 
increase  the  fidelity  of  the  results  significantly.  A  mesh  with  large  numbers  of  elements 
would  require  longer  computational  time  as  well. 

To  achieve  smaller  Y-plus  values,  engineers  should  refine  the  mesh  parameters  at 
critical  area  of  interest.  This  will  improve  the  fidelity  of  the  simulation  without  having  to 
increase  the  mesh  size  significantly.  For  this  report,  the  areas  of  interest  are  the  surfaces 
of  the  projectile  especially  at  the  nose  tip  and  base.  Inflation  layers  can  be  set  at  the 
surface  of  the  projectile  for  refinement  as  shown  in  Figure  7.  The  maximum  thickness  for 
the  inflation  layer  will  depend  on  the  projectile  geometry  size  and  the  number  of  layers 
that  can  be  inserted  will  depend  on  the  maximum  thickness  and  growth  rate  defined. 
Specific  face  sizing  can  be  applied  at  surfaces  of  the  base  and  nose  tip  to  further  refine 
the  mesh  as  shown  in  Figure  8  and  9. 


Figure  7.  Inflation  Layers  (CFX). 
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Figure  8.  Face  Sizing  at  Base  (CFX). 


Figure  9.  Face  Sizing  at  Nose  Tip  (CFX). 


3.  Pre-processor 

The  model  was  set  to  analysis  the  flow  at  steady-state  conditions  with  the  fluid 
definition  set  as  air-ideal  gas.  The  simulations  were  run  with  both  the  k-  s  and  SST  model 
for  high  speed  conditions  and  compressibility  effects.  For  speed  of  MO. 5  to  Ml.O  and 
Ml. 5  to  M2. 6,  the  inlet  flow  conditions  were  set  as  mixed  flow  and  supersonic  flow 
respectively  with  a  static  temperature  setting  of  288. 15K.  Outlet  conditions  were  set  as 
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subsonic  flow  throughout  with  a  static  pressure  of  0  Pa.  Openings  were  set  as 
entrainments  with  a  static  pressure  of  0  Pa  as  well. 

4.  Post-processor 

Flow  across  the  projectile  was  studied  by  inserting  a  plane  across  the  X-Y  axis  on 
the  symmetry  plane.  Velocity  vectors  with  local  Mach  number  and  pressure  contours  as 
shown  in  Figure  10,  1 1  and  12  can  be  applied  onto  the  plane  to  understand  the  flow  field 
of  the  projectile  in  flight. 


Figure  10.  Velocity  Vector  (CFX). 


The  function  calculator  was  used  to  extract  the  respective  Y-plus,  force  and 
torque  results  from  the  simulations  as  shown  in  Appendix  D. 

C.  AP09 

1.  Introduction 

AP09  uses  a  semi-empirical  technique  for  the  estimation  of  aerodynamic 
properties  on  spin  stabilized  projectiles  as  well  as  various  missile  configurations.  The 
program  was  developed  in  1972  and  has  since  evolved  with  the  combination  of 
experimental  results  with  theoretical  methodology  to  produce  reasonable  aerodynamic 
prediction  for  flight  bodies. 

In  this  study,  the  AP09  was  used  to  compute  the  various  aerodynamic  properties 
of  the  M3 3  projectile  to  be  compared  with  the  results  derived  from  CFX.  Details  of  the 
processes  and  input  parameters  can  be  found  in  Appendix  C. 

2.  M33  Projectile  Configuration  and  Input  Parameters 

Table  1  shows  the  input  parameters  for  the  M3 3  projectile  configuration  with 
reference  to  [14]. 
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General 

Reference  Diameter 

12.95mm 

Distance  of  Moment  from  Nose 

Tip 

33.15mm 

Nose  Geometry 

Nose  Cross  Sectional  Shape 

Circular 

Circular  Radius  at  Nose 

6.475mm 

Nose  Profile 

Secant  Ogive  Blunt 

Radius  of  Curvature 

113.57mm 

Spherical  Cap  Radius 

2.33mm 

Length  of  Nose 

33.15mm 

Afterbody  Geometry 

Longitudinal  Afterbody 

Coordinate  from  Nose  Tip 

47.66mm 

Boattail/  Flare 

Longitudinal  Boattail/  Flare 
Coordinate  from  Nose  Tip 

57.76mm 

Corresponding  Boattail/  Flare 
Characteristic  Half  Width 

4.875mm 

Options 

Sea  level 

Dynamic  derivatives 

Spin  stabilized 

Table  1.  M33  Input  Parameters  (AP09). 


3.  Results 

The  program  was  run  for  “Mach  Sweep”  from  M=  0.5  to  3.0  with  an  interval  size 
of  0.15  at  0°  and  2°  angle  of  attack.  The  results  obtained  were  compared  to  previous 
results  in  the  API  reference  [3]  to  check  for  error  rate.  The  results  obtained  agree  fairly 
well  with  the  API  reference  as  shown  in  Figure  13  and  14. 
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Figure  13.  Cdo  vs  Mach  (AP09). 


Figure  14.  Cno  vs  Mach  (AP09). 
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From  the  same  configuration  and  input  parameters,  the  plots  for  Cl  vs  a  and  Cd 
vs  a  at  M=  0.5  for  0°  to  85°  angles  of  attack  as  shown  in  Figure  15  and  16  were  generated 
for  comparison  with  the  results  from  CFX. 


20 


V.  RESULTS  AND  ANALYSIS 


A.  LINEARITY  AT  LOW  ANGLES  OE  ATTACK 

To  calculate  CNa  and  Cmo,  we  must  first  ensure  that  the  aerodynamic  properties 
are  linear  at  low  angles  of  attack.  Figure  18  and  20  shows  that  the  results  for  Cn  are  fairly 
linear  for  subsonie  and  supersonic  Mach  numbers  up  to  6°  angles  of  attack.  For  transonic 
Mach  numbers  as  shown  in  Figure  19,  Cn  is  only  linear  up  to  4°  angles  of  attack. 
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B. 


EFFECTS  OF  MACH  NUMBERS 


1.  k-  epsilon  Model 

The  results  from  the  k-  s  model  show  good  agreement  with  experimental  data  for 
both  Cno  and  Cmu,  but  tend  to  over  predict  for  Cdo  at  higher  Mach  numbers  as  shown  in 
Figures  21,  22  and  23.  This  could  be  due  to  the  models  insufficiency  in  cases  with  large 
adverse  pressure  gradient  and  separated  flow,  as  stated  in  [12], 
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Figure  22. 


2.  Shear  Stress  Transport  Model 

The  SST  model  on  the  other  hand  produces  good  agreement  with  experimental 
data  only  for  supersonic  Mach  numbers  as  shown  in  Figures  24,  25  and  26.  Although  it 
tends  to  over  predict  Cdo  as  well,  the  results  for  M=  1.5  to  2.6  are  closer  to  the 
experimental  data  as  compared  to  the  k-  s  model.  The  poor  results  for  the  subsonic  and 
transonic  Mach  numbers  especially  at  M=1.0  could  be  due  to  the  models  tendency  to 
over-predict  the  turbulence  level  in  regions  with  large  normal  strain  and  acceleration  as 
stated  in  [12]. 
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Figure  23.  Cno  vs  Mach  Numbers_SST  Model  (CFX). 
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Figure  24.  CMa  vs  Mach  Numbers_SST  Model  (CFX). 
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3. 


Combined  k-  £  and  SST  Model 


Combining  the  results  from  the  k-  s  model  for  M=  0.5  to  1.0  and  STT  model  for 
M=1.5  to  2.6  produees  good  agreement  against  the  experimental  range  data  for  CNa  and 
Cmo  as  shown  in  Figures  27  and  28.  Cdo  tends  to  be  over  predieted  by  an  average  error  of 
11.6%  with  the  highest  error  oceurring  at  M=  1.5  as  shown  in  Figure  29.  Detailed  data 
and  flow  diagrams  of  the  projectile  at  various  Mach  numbers  can  be  found  in  Appendix 
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Figure  26.  CNa  vs  Mach  Numbers_Combined  k-  s/  SST  Model  (CFX). 
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c. 


EFFECTS  OF  ANGLES  OF  ATTACK 


The  model  was  ran  for  various  angles  of  attaek  at  M=  0.5  to  prediet  the  stall  value 
for  the  M3  3  projeetile.  The  results  from  the  CFX  eode  were  eompared  with  the  results 
obtained  from  AP09  as  shown  in  Figures  30  and  3 1 .  Although  both  eodes  predieted  the 
stall  value  to  occur  at  60°  angle  of  attack,  the  results  differ  in  magnitudes  and  agree  only 
up  to  20°.  Detailed  data  and  flow  diagrams  of  the  projectile  at  various  angles  of  attack 
can  be  found  in  Appendix  D. 
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The  differences  in  magnitude  sparked  off  an  investigation  to  study  the  flow  over 
the  cross-section  of  the  projectile  body.  The  model  was  run  at  M=  0.5  for  45°,  75°  and 
85°  angle  of  attack  with  the  flow  diagrams  shown  in  Appendix  D.  From  Figure  31,  it  was 
understood  that  a  limitation  to  the  half-bodied  symmetry  model  was  its  insufficiency  to 
account  for  the  effects  across  the  Y-Z  plane.  You  could  see  that  the  flows  over  the 
symmetry  plane  were  forced  to  re-circulate  back  as  if  it  has  encountered  a  wall.  This 
phenomenon  is  unnatural  as  flows  may  not  necessary  be  symmetrical  and  could  even  be 
shed  away.  Errors  due  to  this  limitation  may  be  magnified  at  higher  angles  of  attack  when 
the  effects  of  the  cross-sectional  flow  also  occur  at  the  nose. 
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Figure  3 1 .  Cross  Sectional  Flow  (CFX) 


D.  SENSITIVITY  STUDIES  FOR  CFX 

1.  Flow  vs  Orientation  Adjustment  Approach 

The  simulation  was  run  using  a  half-bodied  symmetry  model  and  by  adjusting  the 
flow  component  at  various  angles  of  attack.  A  sensitivity  study  on  adjusting  the 
orientation  of  the  projectile  instead  of  the  flow  was  done  to  compare  the  results.  The  body 
of  the  projectile  was  tilted  by  2°  and  run  with  the  k-  s  model.  The  results  from  the 
simulation  as  shown  in  Figure  32  shows  fairly  good  agreement  with  the  results  obtained 
from  the  flow  adjustment  approach. 
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2.  Control  Volume 

In  actual  wind  tunnel  tests,  measures  need  to  be  taken  to  minimize  the  effeets  of 
re-circulated  flow  from  the  walls  to  the  projectile.  This  measure  was  taken  into  aeeount  in 
the  CFX  simulation  by  using  a  control  volume  that  is  10  times  the  projeetile  length.  An 
additional  sensitivity  study  was  earned  out  to  verify  that  this  eontrol  volume  size  was 
suffieient  for  proper  simulation.  A  control  volume  of  20  times  the  projeetile  length  was 
used  with  the  SST  model  and  the  results  agree  fairly  well  with  the  previous  eontrol 
volume  as  shown  in  the  Figures  33  and  34. 
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VI.  CONCLUSION  AND  RECOMMENDATIONS 


The  CFX  code,  based  on  the  numerical  solution  of  the  full  Navier-Stokes 
equations,  took  about  six  hours  to  compute  for  one  data  point  using  a  mesh  size  of  over 
650,000  nodes  and  2  million  elements.  Mesh  refinement  using  inflation  layers  and  face 
sizing  at  critical  areas  produces  better  simulation  results  and  higher  fidelity. 

The  combined  k-  s  and  SST  turbulence  model  produces  reasonably  good  results 
that  showed  good  agreement  with  experimental  range  data  for  Cnu  and  Cmu-  Cdo  tends  to 
be  over  predicted  by  an  average  error  of  11.6%  with  the  highest  error  occurring  at  M= 
1.5.  As  the  simulation  was  run  using  the  free  flight  motion  of  the  M33  projectile,  a  large 
source  of  error  may  be  contributed  due  to  the  rotational  effects  of  the  spin-stabilized 
round  during  range  testing.  For  simulation,  it  is  recommended  to  use  the  k-  s  model  for 
subsonic  Mach  numbers  that  are  M<  1.0  and  the  SST  model  for  supersonic  Mach 
numbers  that  are  M>  1.0. 

Stall  value  of  the  M33  projectile  was  predicted  at  60°  angle  of  attack.  Limitations 
to  the  cross  flow  analysis  of  the  half-bodied  symmetry  model  require  further  studies  on  a 
full-bodied  model  for  better  fidelity. 

Sensitivity  studies  showed  that  a  control  volume  of  10  times  the  projectile  length 
is  sufficient  to  avoid  errors  due  to  re-circulation  of  flows  from  the  walls.  Computational 
approach  of  adjusting  the  flow  at  various  angles  of  attack  also  shows  good  agreement 
with  the  results  from  adjusting  the  projectile  orientation  for  various  angles  of  attack. 

In  conclusion,  this  report  serves  as  a  further  study  from  previous  report  [4]  in 
using  the  CFX  code  for  the  prediction  of  aerodynamic  properties  of  projectiles  at  various 
angles  of  attack.  This  report  also  serves  as  a  foundation  for  further  studies  into  the  cross 
flow  analysis  and  rotational  effects  of  the  spin-stabilized  M3 3  projectile. 
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APPENDIX  A.  SOLIDWORKS  PROCESSES  AND  PROCEDURES 


A.  INTRODUCTION 

Gouge  material  for  SolidWorks  processes  and  procedures  can  be  taken  from  the 
laboratory  work  in  [  1 5] . 

B.  2D  MODEL 

Select  the  front  plane  (XY  plane)  as  the  reference  plane.  Plot  the  individual  points 
of  the  projectile  on  the  XY  axis  at  coordinates  (0,0),  (0,  4.875E-3),  (lO.lOlE-3,  6.475E- 
3),  (24.605E-3,  6.475E-3),  (55.426E-3,  2.33 lE-3)  and  (57.757E-3,  0). 


Figure  35.  Plotting  Individual  Points  (SolidWorks). 


Connect  the  individual  points  using  straight  lines  for  the  afterbody  and  tangent 
arcs  for  the  nose  cone  and  nose  tip.  Set  the  radius  for  the  nose  cone  and  nose  tip  to  be 
113.5715E-3  and2.331E-3  respectively. 
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Figure  36.  Straight  Lines  for  Afterbody  (SolidWorks). 


Figure  37.  Tangent  Arcs  for  nose  (SolidWorks). 


The  control  volume  that  is  10  times  the  size  of  the  reference  projectile  is  drawn 
using  a  rectangular  box  on  the  XY  axis  with  coordinates  (-570E-3,  570E-3),  (-570E-3,  - 
570E-3),  (600E-3,  570E-3)  and  (600E-3,  -570E-3). 
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C.  3D  MODEL 

The  3D  model  of  the  control  volume  is  generated  by  extruding  the  2D  control 
volume  in  the  Z  axis  by  -570E-3. 


Figure  39.  3D  Control  Volume  (SolidWorks). 

Perform  a  revolved  cut  for  360°  with  the  2D  sketch  of  the  projectile  and  it  will 

leave  a  cavity  of  the  half-bodied  projectile  in  the  control  volume  for  simulation  purposes. 
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Figure  40.  Revolved  Cut  for  2D  Projectile  Sketch  (SolidWorks). 
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APPENDIX  B.  CFX  PROCESSES  AND  PROCEDURES 


A.  INTRODUCTION 

Gouge  material  for  CFX  processes  and  procedures  can  be  taken  from  the 
laboratory  work  in  [  1 5] . 

B.  GEOMETRY 

After  starting  up  the  ANYSYS-CFX,  import  the  M3  3  geometry  fde  from 
SolidWorks  that  was  saved  as  a  parasolid  fde. 


3 

^  Mesh 

(Sii) 

New  Geometry... 

4 

Setup 

Import  Geometry  ► 

^"1  Browse... 

5 

^  Solution 

Duplicate 

rS  M33.x_t 

6 

1^  Results 

Transfer  Data  From  New  ► 

cheeheng.SLDPRT 

Fluid  Flow  (CFX; 

Transfer  Data  To  New  ► 

airfoil.SLDPRT 

Update 

qS.SLDPRT 

Reset 

Rename 

Properties 

Quick  Flelp 

Figure  4 1 .  Geometry  Interface  (CFX). 


C.  MESH 

I.  Named  Selections 

Inside  the  mesh  interface,  insert  and  define  named  selections  under  geometry. 
Define  the  projectile  surfaces  as  “body”  and  the  rest  of  the  control  volume  surfaces  as 
shown  in  Figure  42. 
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Figure  42.  Named  Selections  (CFX). 


2.  Mesh  Sizing 


Table  2  shows  the  input  parameters  that  were  set  for  the  mesh. 


Sizing 

Relevance  Center 

Fine 

Curvature  Angle 

2° 

Min  Size 

1  X  10'"^m 

Max  Size 

0.1m 

Max  Face  Size 

0.1m 

Inflation 

Use  Automatic  Inflation 

All  faces  in  chosen  selection 

Named  Selection 

Body 

Inflation  Option 

Total  Thickness 

Number  of  Layers 

12 

Growth  Rate 

0.8 

Max  Thickness 

1  X  10‘^m 

Table  2.  Mesh  Sizing  Input  Parameters  (CFX). 


3.  Face  Sizing 

Insert  face  sizing  under  mesh  and  select  the  base  and  nose  tip  surface.  Select  type 
to  be  element  size  of  1x10'^  m. 
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Figure  43.  Face  Sizing  (CFX). 


Insert  method  under  mesh  to  be  automatic  and  the  program  should  generate  a 
mesh  with  over  650,000  nodes  and  2  million  elements. 

D.  PRE-PROCESSOR 

1.  Domain  Conditions 

Set  the  analysis  type  to  steady-state.  Inside  the  default  domain,  define  material  for 
fluid  and  particle  under  basic  setting  to  be  air-ideal  gas.  Set  the  fluid  model  for  heat 
transfer  to  be  for  total  energy  and  the  turbulence  model  to  be  either  SST  of  k-  s  model. 
For  SST  model,  insert  additional  parameter  under  transitional  turbulence  for  Gamma 
Theta  Model.  Check  the  option  for  high  speed  (compressible)  wall  heat  transfer  model. 
Initialize  the  default  domain  condition  as  the  Ifee-stream  condition  with  Cartesian 
velocity  values  set  for  the  u  and  v  component  depending  on  the  Mach  number  and  angle 
of  attack. 

2.  Boundary  Conditions 

Define  and  insert  boundary  conditions  under  the  default  domain  as  the  following 


parameters. 


Boundary 

Basic  Setting 

Boundary  Details 

body 

Boundary  Type  -  Wall 
Location  -  body 

Mass  and  Momentum  -  No  Slip  Wall 

Wall  Roughness -Smooth  Wall 

Heat  Transfer  -  Adiabatic 
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Boundary 

Basic  Setting 

Boundary  Details 

inlet 

Boundary  Type  -  Inlet 
Location  -  inlet 

Flow  Regime  -  Mixed  Flow  (MO. 5  to  Ml.O)/ 
Supersonic  Flow  (Ml. 5  to  M2. 6) 

Turbulence  -  Medium  (Intensity  =  5%) 

Fleat  Transfer  -  Static  Temperature  (288. 15K) 

outlet 

Boundary  Type  -  Outlet 
Location  -  outlet 

Flow  Regime  -  Subsonic 

Mass  and  Momentum  -  Static  Pressure  (OPa) 

top 

Boundary  Type  -  Opening 
(a=0°)/  Outlet  (a=2°) 

Location  -  top 

For  Opening: 

Flow  Regime  -  Subsonic 

Mass  and  Momentum  -  Entrainment  (OPa) 

Turbulence  -  Zero  Gradient 

Heat  Transfer  -  Static  Temperature  (288. 15K) 

For  Outlet:  Same  as  outlet 

bottom 

Boundary  Type  -  Opening 
(a=0°)/  Inlet  (a=2°) 

Location  -  bottom 

For  Opening: 

Flow  Regime  -  Subsonic 

Mass  and  Momentum  -  Entrainment  (OPa) 

Turbulence  -  Zero  Gradient 

Heat  Transfer  -  Static  Temperature  (288. 15K) 

For  Inlet:  Same  as  inlet 

side 

Boundary  Type  -  Opening 
Location  -  side 

Flow  Regime  -  Subsonic 

Mass  and  Momentum  -  Entrainment  (OPa) 

Turbulence  -  Zero  Gradient 

Heat  Transfer  -  Static  Temperature  (288. 15K) 

sym 

Boundary  Type  -  Symmetry 
Location  -  side 

Table  3.  Boundary  Condition  Input  Parameters  (CFX). 
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Figure  44.  Boundary  Conditions  at  0°  (CFX). 


3.  Solver 

Insert  expert  parameters  under  solver  and  set  eonvergenee  eontrol  in  high  speed 
models  for  three  eontinuity  loops.  For  solver  control,  set  the  turbulence  numeric  to  high 
resolution.  Set  maximum  number  of  iterations  to  be  100  with  a  residual  target  of  1x10'^. 
In  the  advance  options,  set  high  speed  numeric  for  compressibility  control. 

E.  POST-PROCESSOR 

1.  Flow  Field 

Insert  Plane  on  the  XY-axis  and  set  the  Z-axis  distance  to  be  -0.05x10'^  m.  Insert 
vector  on  plane  with  variable  set  as  velocity  to  obtain  the  velocity  field  vector.  Insert 
contour  on  plane  and  set  variable  as  local  Mach  number  or  pressure  to  obtain  the  Mach 
number  and  pressure  contours. 

2.  Forces 

Use  the  function  calculator  to  obtain  the  values  for  the  forces  acting  on  the 

projectile.  Table  4  shows  the  input  parameters  to  obtain  the  results. 
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Y-Plus 

Function  -  Ave 

Location  -  Default  Domain 

Variable  -  Y-plus 

Forces 

Function  -  Force 

Location  -  Body 

Direction  -  X/  Y/  Z 

Torque 

Function  -  Torque 

Location  -  Body 

Axis  -  X/  Y/  Z 

Table  4.  Post-Processor  Input  Parameters  (CFX). 
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APPENDIX  C.  AP09  PROCESSES  AND  PROCEDURES 


A.  M33  PROJECTILE  CONFIGURATION  AND  INPUT  PARAMETERS 

To  input  the  geometry  and  aerodynamic  parameters  for  the  M3  3,  create  a  new 
“body-alone”  configuration  and  configure  the  geometry  input  units  to  millimeters. 


Figure  45.  Configuration-New-Body- Alone  (AP09). 


Figure  46.  Inputs-Geometry-Geometry  (AP09). 
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B.  GEOMETRY  INPUTS 

Inside  the  geometry  data  entry  interface,  set  the  following  body  roughness  for 
typical  flight  configuration,  reference  diameter  or  width  of  the  body  to  be  12.95mm  and 
distance  of  moment  reference  from  nose  tip  to  be  33.15mm. 


Figure  47.  Geometry  Data  Entry  Interface  (AP09). 


To  configure  nose  profile,  select  “Nose  Geometry”  under  the  geometry  data  entry 
interface.  Select  nose  cross  sectional  shape  to  be  circular  and  set  circular  radius  to 
6.475mm.  Next,  select  “Secant  Ogive  Blunt”  for  nose  profile  and  set  radius  of  curvature 
to  be  1 13.57mm,  spherical  cap  radius  to  be  2.33mm  and  length  of  nose  tip  to  shoulder  to 
be  33.15mm. 
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Figure  48.  Nose  Geometry  Data  Entry  Interface  (AP09). 


Figure  49. 


IBS 


Nose  Profile  Data  Entry  Interface  (AP09). 


To  configure  the  body  geometry  of  the  M33,  select  “Afterbody  Geometry”  under 
the  geometry  data  entry  interface  and  set  longitudinal  afterbody  coordinate  from  nose  tip 
to  be  47.66mm. 
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Figure  50.  Afterbody  Geometry  Data  Entry  Interface  (AP09). 


^11 1  Airi?  il-  I  : 


Figure  5 1 .  Afterbody  Profde  Data  Entry  Interface  (AP09). 


Lastly  for  geometry  inputs,  check  “yes”  for  “Boattail/  Flare”  under  the  geometry 
data  entry  interface  and  set  longitudinal  boattail/  flare  coordinate  from  nose  tip  to  be 
57.76mm  and  corresponding  boattail/  flare  characteristic  half  width  to  be  4.875mm. 
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Figure  52.  Boattail/  Flare  Data  Entry  Interface  (AP09). 


C.  AERODYNAMIC  INPUTS 

To  enter  aerodynamic  parameters,  select  “Aerodynamics”  under  “Inputs”  and  set 
the  free-stream  conditions  for  Mach  sweep. 


Figure  53.  Inputs-Aerodynamics  (AP09). 
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Figure  54.  Aerodynamics  Data  Entry  Interface-Free-Stream  Conditions  (AP09). 


Inside  the  Mach  sweep  data  entry  interface,  set  initial  Mach  number  to  be  M=  0.5, 
final  Mach  number  to  be  M=  3,  interval  size  to  be  0.15,  constant  angle  of  attack  to  be  0° 


or  2°  and  altitude  to  be  1  feet  (sea  level). 


Figure  55.  Mach  Sweep  Data  Entry  Interface  (AP09). 


After  keying  in  the  aerodynamic  parameters,  select  “Options”  under  the 
aerodynamic  data  entry  interface.  Check  “yes”  for  dynamic  derivatives,  spin  stabilized 
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and  set  aerodynamic  smoother  to  “on”.  A  title  is  required  in  order  to  proceed.  Simply 
enter  CASEl  or  anything  of  preference. 


Figure  56.  Aerodynamics  Data  Entry  Interface-Options  (AP09). 


Figure  57.  Aerodynamic  Options  Data  Entry  Interface  (AP09). 

D.  RESULTS 

The  geometry  sketch  from  the  geometry  inputs  can  be  previewed  under 
“Generate”.  After  confirming  the  geometry,  generate  the  aerodynamic  output  file.  This 
process  will  take  a  few  seconds. 
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Figure  58.  Generate  -  Aerodynamic  Output  File  (AP09). 


Select  “Plots”  under  “Output”  to  retrieve  the  generated  results  for  the  M33 
projectiles.  Tabulated  results  from  the  output  files  can  be  saved  as  text  documents  and 
expressed  using  excel  spreadsheet. 


Figure  59.  Outputs-Aerodynamics-Plots/  Tables  (AP09). 
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Figure  60.  Output  Interface  (AP09). 
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APPENDIX  D.  DETAILED  RESULTS 


A.  CFX 


1.  Simulation  Data 


Simulation  Data  for  Cn  vs  AOA  (MO.5) 

Drag  Ref  Area  (m^) 

0.00013 

Density  (kg/m^) 

1.180 

Mach  Number 

0.500 

Velocity  (m/s) 

170.150 

an 

-4.000 

-2.000 

0.000 

2.000 

4.000 

6.000 

a  (rad) 

-0.070 

-0.035 

0.000 

0.035 

0.070 

0.105 

u  (m/s) 

169.735 

170.046 

170.150 

170.046 

169.735 

169.218 

V  (m/s) 

-11.871 

-5.939 

0.000 

5.939 

11.871 

17.788 

Fy(N) 

-0.107 

-0.060 

-0.001 

0.069 

0.137 

0.198 

F.(N) 

0.259 

0.253 

0.138 

0.155 

0.261 

0.265 

Fl(N) 

-0.089 

-0.051 

-0.001 

0.064 

0.118 

0.169 

Fd  (N) 

0.265 

0.255 

0.138 

0.158 

0.269 

0.285 

Cl 

-0.080 

-0.046 

-0.001 

0.057 

0.106 

0.152 

Cd 

0.239 

0.230 

0.125 

0.142 

0.243 

0.256 

Cn 

-0.097 

-0.054 

-0.001 

0.062 

0.123 

0.178 

Table  5.  Simulation  Results  for  Cn  vs  a  Subsonic  (CFX). 


Simulation  Data  for  Cn  vs  AOA  (M0.8) 

Drag  Ref  Area  (m^) 

0.00013 

Density  (kg/m^) 

1.180 

Mach  Number 

0.800 

Velocity  (m/s) 

272.240 

an 

-4.000 

-2.000 

0.000 

2.000 

4.000 

6.000 

a  (rad) 

-0.070 

-0.035 

0.000 

0.035 

0.070 

0.105 

u  (m/s) 

271.577 

272.074 

272.240 

272.074 

271.577 

270.748 

V  (m/s) 

-18.993 

-9.502 

0.000 

9.502 

18.993 

28.461 

Fy(N) 

-0.177 

-0.093 

0.000 

0.166 

0.340 

0.393 

F.(N) 

0.580 

0.572 

0.469 

0.475 

0.600 

0.624 

Fl(N) 

-0.136 

-0.073 

0.000 

0.150 

0.298 

0.326 
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Simulation  Data  for  Cn  vs  AOA  (MO.8) 

Fd  (N) 

0.591 

0.575 

0.469 

0.481 

0.622 

0.662 

Cl 

-0.048 

-0.026 

0.000 

0.053 

0.105 

0.115 

Cd 

0.208 

0.202 

0.165 

0.169 

0.219 

0.233 

Cn 

-0.062 

-0.033 

0.000 

0.059 

0.120 

0.138 

Table  6.  Simulation  Results  for  Cn  vs  a_Transonie  (CFX). 


Simulation  Data  for  Cn  vs  AOA  (Ml. 5) 

Drag  Ref  Area  (m^) 

0.00013 

Density  (kg/m^) 

1.180 

Mach  Number 

1.500 

Velocity  (m/s) 

510.450 

an 

-4.000 

-2.000 

0.000 

2.000 

4.000 

6.000 

a  (rad) 

-0.070 

-0.035 

0.000 

0.035 

0.070 

0.105 

u  (m/s) 

509.206 

510.139 

510.450 

510.139 

509.206 

507.653 

V  (m/s) 

-35.612 

-17.817 

0.000 

17.817 

35.612 

53.363 

Fy(N) 

-1.247 

-0.626 

0.000 

0.863 

1.756 

2.706 

F.(N) 

4.316 

4.242 

4.194 

4.240 

4.346 

4.468 

Fl(N) 

-0.942 

-0.478 

0.000 

0.715 

1.449 

2.224 

Fd  (N) 

4.393 

4.262 

4.194 

4.267 

4.458 

4.727 

Cl 

-0.094 

-0.048 

0.000 

0.072 

0.145 

0.223 

Cd 

0.440 

0.426 

0.420 

0.427 

0.446 

0.473 

Cn 

-0.125 

-0.063 

0.000 

0.086 

0.176 

0.271 

Table  7.  Simulation  Results  for  Cn  vs  a_Supersonic  (CFX). 


Simulation  Data  (k-e  Model) 

Diameter 

(m) 

0.01295 

Drag  Ref 
Area  (m^) 

0.00013 

Air  Density 
(kg/m’) 

1.180 

length  (m) 

0.058 

XcG  from 
Base  (m) 

0.023 

Mach 

Number 

0.500 

0.800 

1.000 

1.500 

1.950 

2.600 

Velocity 

(m/s) 

170.150 

272.240 

340.300 

510.450 

663.585 

884.780 

an 

0.000  2.000 

0.000  2.000 

0.000  2.000 

0.000  2.000 

0.000  2.000 

0.000  2.000 
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Simulation  Data  (k-e  Model) 

a  (rad) 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

u  (m/s) 

170.1 

50 

170.0 

46 

272.2 

40 

272.0 

74 

340.3 

00 

340.0 

93 

510.4 

50 

510.1 

39 

663.5 

85 

663.1 

81 

884.7 

80 

884.2 

41 

V  (m/s) 

0.000 

5.939 

0.000 

9.502 

0.000 

11.87 

8 

0.000 

17.81 

7 

0.000 

23.16 

2 

0.000 

30.88 

2 

Y+ 

0.926 

1.650 

2.030 

2.800 

3.250 

3.700 

Fy(N) 

-0.001 

0.069 

0.000 

0.166 

-0.001 

0.290 

0.001 

0.865 

0.005 

1.597 

-0.005 

2.893 

Fx(N) 

0.138 

0.155 

0.469 

0.475 

1.534 

1.549 

4.228 

4.285 

6.287 

6.341 

9.168 

9.212 

Fl  (N) 

-0.001 

0.064 

0.000 

0.150 

-0.001 

0.235 

0.001 

0.715 

0.005 

1.374 

-0.005 

2.569 

Fd  (N) 

0.138 

0.158 

0.469 

0.481 

1.534 

1.558 

4.228 

4.313 

6.287 

6.393 

9.168 

9.307 

Tz  (Nm) 

0.000 

0.003 

0.000 

0.010 

0.000 

0.015 

0.000 

0.031 

0.000 

0.049 

0.000 

0.075 

Cl 

-0.001 

0.057 

0.000 

0.052 

0.000 

0.052 

0.000 

0.071 

0.000 

0.080 

0.000 

0.084 

Cd 

0.123 

0.140 

0.163 

0.167 

0.341 

0.346 

0.418 

0.426 

0.367 

0.374 

0.301 

0.306 

Cn 

-0.001 

0.062 

0.000 

0.058 

0.000 

0.064 

0.000 

0.085 

0.000 

0.093 

0.000 

0.095 

Cm 

-0.001 

0.127 

0.000 

0.147 

0.000 

0.148 

0.000 

0.131 

0.000 

0.124 

0.000 

0.107 

Cnu 

1.794 

1.660 

1.852 

2.445 

2.664 

2.729 

Cmu 

3.646 

4.217 

4.241 

3.762 

3.539 

3.053 

Table  8.  Simulation  Results  for  k-  s  model  (CFX). 


Simulation  Data  (SST  Model) 

Diameter 

(m) 

0.01295 

Drag  Ref 
Area  (m^) 

0.00013 

Air  Density 
(kg/m’) 

1.180 

length  (m) 

0.058 

XcG  from 
Base  (m) 

0.023 

Mach 

Number 

0.500 

0.800 

1.000 

1.500 

1.950 

2.600 

Velocity 

(m/s) 

170.150 

272.240 

340.300 

510.450 

663.585 

884.780 

an 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

a  (rad) 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

u  (m/s) 

170.1 

50 

170.0 

46 

272.2 

40 

272.0 

74 

340.3 

00 

340.0 

93 

510.4 

50 

510.1 

39 

663.5 

85 

663.1 

81 

884.7 

80 

884.2 

41 

V  (m/s) 

0.000 

5.939 

0.000 

9.502 

0.000 

11.87 

8 

0.000 

17.81 

7 

0.000 

23.16 

2 

0.000 

30.88 

2 

Y+ 

1.500 

2.030 

2.430 

2.800 

3.050 

3.280 

Fy(N) 

-0.002 

0.076 

-0.004 

0.180 

0.000 

0.406 

0.000 

0.863 

0.006 

1.593 

-0.007 

2.872 

Fx(N) 

0.249 

0.257 

0.566 

0.583 

2.430 

2.453 

4.194 

4.240 

6.119 

6.173 

8.790 

8.838 

Fl  (N) 

-0.002 

0.067 

-0.004 

0.160 

0.000 

0.320 

0.000 

0.715 

0.006 

1.376 

-0.007 

2.562 
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Simulation  Data  (SST  Model) 

Fd  (N) 

0.249 

0.259 

0.566 

0.589 

2.430 

2.466 

4.194 

4.267 

6.119 

6.224 

8.790 

8.933 

Tz  (Nm) 

0.000 

0.003 

0.000 

0.010 

0.000 

0.018 

0.000 

0.030 

0.000 

0.047 

0.000 

0.073 

Cl 

-0.002 

0.060 

-0.001 

0.056 

0.000 

0.071 

0.000 

0.071 

0.000 

0.080 

0.000 

0.084 

Cd 

0.221 

0.230 

0.197 

0.204 

0.540 

0.548 

0.414 

0.421 

0.358 

0.364 

0.289 

0.294 

Cn 

-0.002 

0.068 

-0.001 

0.063 

0.000 

0.090 

0.000 

0.085 

0.000 

0.093 

0.000 

0.094 

Cm 

-0.002 

0.132 

-0.001 

0.156 

0.000 

0.169 

0.000 

0.129 

0.000 

0.119 

0.000 

0.104 

Cnu 

1.997 

1.835 

2.585 

2.443 

2.656 

2.711 

Cmu 

3.844 

4.493 

4.850 

3.696 

3.414 

2.988 

Table  9.  Simulation  Results  for  SST  Model  (CFX). 


Simulation  Data  for  Various  a 

Area  (m^) 

0.00013 

Density 

(kg/m’) 

1.180 

Mach 

Number 

0.500 

Velocity 

(m/s) 

170.150 

«(°) 

-4.000 

-2.000 

0.000 

2.000 

4.000 

6.000 

10.000 

15.000 

20.000 

25.000 

30.000 

a  (rad) 

-0.070 

-0.035 

0.000 

0.035 

0.070 

0.105 

0.175 

0.262 

0.349 

0.436 

0.524 

u  (m/s) 

169.73 

5 

170.04 

6 

170.15 

0 

170.04 

6 

169.73 

5 

169.21 

8 

167.56 

4 

164.35 

1 

159.88 

6 

154.20 

4 

147.34 

8 

V  (m/s) 

11.871 

-5.939 

0.000 

5.939 

11.871 

17.788 

29.550 

44.044 

58.202 

71.917 

85.085 

Fy(N) 

-0.107 

-0.060 

-0.001 

0.069 

0.137 

0.198 

0.337 

0.536 

0.742 

0.961 

1.185 

F.(N) 

0.259 

0.253 

0.138 

0.155 

0.261 

0.265 

0.275 

0.293 

0.300 

0.279 

0.238 

Fl  (N) 

-0.089 

-0.051 

-0.001 

0.064 

0.118 

0.169 

0.284 

0.442 

0.595 

0.753 

0.908 

Fd  (N) 

0.265 

0.255 

0.138 

0.158 

0.269 

0.285 

0.329 

0.422 

0.535 

0.659 

0.799 

Cl 

-0.080 

-0.046 

-0.001 

0.057 

0.106 

0.152 

0.255 

0.398 

0.536 

0.678 

0.817 

Cd 

0.239 

0.230 

0.125 

0.142 

0.243 

0.256 

0.297 

0.380 

0.482 

0.594 

0.720 

a(') 

35.000 

40.000 

45.000 

50.000 

55.000 

60.000 

65.000 

70.000 

75.000 

80.000 

85.000 

a  (rad) 

0.611 

0.698 

0.786 

0.873 

0.960 

1.047 

1.135 

1.222 

1.309 

1.396 

1.484 

u  (m/s) 

139.37 

1 

130.33 

3 

120.30 

2 

109.35 

6 

97.577 

85.055 

71.886 

58.169 

44.010 

29.516 

14.797 

V  (m/s) 

97.605 

109.38 

2 

120.32 

6 

130.35 

5 

139.39 

1 

147.36 

6 

154.21 

9 

159.89 

8 

164.36 

0 

167.57 

0 

169.50 

5 

Fv(N) 

1.410 

1.636 

1.864 

2.091 

2.307 

2.500 

2.644 

2.709 

2.661 

2.513 

2.419 

F«(N) 

0.173 

0.086 

-0.017 

-0.137 

-0.264 

-0.392 

-0.508 

-0.594 

-0.572 

-0.419 

-0.339 
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Simulation  Data  for  Various  a 

Fl(N) 

1.056 

1.197 

1.330 

1.449 

1.540 

1.590 

1.577 

1.485 

1.240 

0.849 

0.548 

Fd(N) 

0.951 

1.118 

1.306 

1.514 

1.738 

1.969 

2.181 

2.342 

2.422 

2.402 

2.381 

Cl 

0.951 

1.078 

1.198 

1.305 

1.387 

1.432 

1.421 

1.337 

1.117 

0.764 

0.493 

Cd 

0.856 

1.007 

1.176 

1.364 

1.566 

1.774 

1.965 

2.110 

2.182 

2.163 

2.144 

Table  10.  Simulation  Results  for  Cl  and  Cd  at  High  AOA  (CFX). 


Sensitivity  Study  on  Tilted  Model  for  2°  AOA 

Diameter  (m) 

0.013 

Drag  Ref  Area  (m^) 

0.000 

Air  Density  (kg/m^) 

1.180 

Mach  Number 

0.500 

1.500 

2.600 

Velocity  (m/s) 

170.150 

510.450 

884.780 

an 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

a  (rad) 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

u  (m/s) 

170.150 

170.046 

510.450 

510.139 

884.780 

884.241 

V  (m/s) 

0.000 

5.939 

0.000 

17.817 

0.000 

30.882 

Fy(N) 

-0.001 

0.069 

0.000 

0.732 

-0.007 

2.695 

Fx(N) 

0.138 

0.155 

4.194 

4.271 

8.790 

8.599 

Cl 

-0.001 

0.061 

0.000 

0.072 

0.000 

0.089 

Cd 

0.123 

0.138 

0.414 

0.422 

0.289 

0.283 

Cn 

-0.001 

0.066 

0.000 

0.087 

0.000 

0.098 

Cnu 

1.914 

2.492 

2.825 

Table  1 1 .  Sensitivity  Study  on  Tilted  Body  vs  Tilted  Flow  (CFX). 


Sensitivity  Study  on  20X  Control  Volume  (SST) 

Diameter  (m) 

0.013 

Drag  Ref  Area  (m^) 

0.000 

Air  Density  (kg/m^) 

1.180 

Mach  Number 

0.500 

0.800 

1.500 

Velocity  (m/s) 

170.150 

272.240 

510.450 

of) 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

a  (rad) 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

u  (m/s) 

170.150 

170.046 

272.240 

272.074 

510.450 

510.139 

V  (m/s) 

0.000 

5.939 

0.000 

9.502 

0.000 

17.817 

Fy(N) 

0.001 

0.075 

0.002 

0.185 

0.000 

0.867 
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Sensitivity  Study  on  20X  Control  Volume  (SST) 

Fx(N) 

0.243 

0.250 

0.552 

0.568 

4.173 

4.222 

Fl  (N) 

0.001 

0.066 

0.002 

0.165 

0.000 

0.719 

Fd  (N) 

0.243 

0.252 

0.552 

0.574 

4.173 

4.250 

Cl 

0.001 

0.059 

0.001 

0.057 

0.000 

0.071 

Cd 

0.216 

0.224 

0.192 

0.199 

0.412 

0.420 

Cn 

0.001 

0.067 

0.001 

0.064 

0.000 

0.086 

Cnu 

1.885 

1.824 

2.451 

Table  12.  Sensitivity  Study  on  20X  Control  Volume  (CFX). 
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2.  Flow  Diagrams  at  Various  Mach  Numbers 
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Velocity  Vector 

Local  Mach  Contour 

Local  Pressure  Contour 

M0.8_0° 
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Velocity  Vector 

Local  Mach  Contour 

Local  Pressure  Contour 

M1.0_0° 
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Velocity  Vector 

Local  Mach  Contour 

Local  Pressure  Contour 

M2.6_0° 

I, 


Figure  61.  Flow  Diagrams  at  Various  Mach  Numbers  for  0°  and  2°  (CFX). 
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3 


Flow  Diagrams  at  Various  Angles  of  Attack 
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Velocity  Vector 

Local  Mach  Contour 

Local  Pressure  Contour 
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Velocity  Vector 

Local  Mach  Contour 

Local  Pressure  Contour 

M0.5_85° 

Figure  62.  Flow  Diagram  at  Various  Angles  of  Attack  (CFX). 
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Cross-Sectional  Flow 
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Figure  63.  Cross-Sectional  Flow  across  Symmetry  Plane  (CFX). 
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B.  AP09 


Simulation  Data  (AP09) 

Mach 

Number 

0.5 

0.65 

0.8 

0.95 

1.1 

1.25 

an 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

a  (rad) 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

Cl 

0.000 

0.050 

0.000 

0.047 

0.000 

0.044 

0.000 

0.042 

0.000 

0.040 

0.000 

0.055 

Cd 

0.097 

0.111 

0.098 

0.112 

0.116 

0.129 

0.194 

0.208 

0.337 

0.358 

0.350 

0.374 

Cn 

0.000 

0.054 

0.000 

0.051 

0.000 

0.048 

0.000 

0.049 

0.000 

0.052 

0.000 

0.068 

Cno 

1.542 

1.457 

1.389 

1.410 

1.503 

1.948 

Mach 

Number 

1.4 

1.55 

1.7 

1.85 

2 

2.15 

an 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

a  (rad) 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

Cl 

0.000 

0.063 

0.000 

0.070 

0.000 

0.076 

0.000 

0.072 

0.000 

0.073 

0.000 

0.073 

Cd 

0.345 

0.369 

0.336 

0.359 

0.330 

0.352 

0.311 

0.331 

0.310 

0.329 

0.309 

0.326 

Cn 

0.000 

0.076 

0.000 

0.082 

0.000 

0.088 

0.000 

0.084 

0.000 

0.084 

0.000 

0.084 

Cno 

2.172 

2.363 

2.528 

2.392 

2.419 

2.416 

Mach 

Number 

2.3 

2.45 

2.6 

2.75 

2.9 

3 

an 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

0.000 

2.000 

a  (rad) 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

0.000 

0.035 

Cl 

0.000 

0.073 

0.000 

0.072 

0.000 

0.071 

0.000 

0.070 

0.000 

0.069 

0.000 

0.068 

Cd 

0.306 

0.321 

0.302 

0.316 

0.299 

0.311 

0.296 

0.307 

0.292 

0.302 

0.289 

0.299 

Cn 

0.000 

0.084 

0.000 

0.083 

0.000 

0.082 

0.000 

0.081 

0.000 

0.079 

0.000 

0.078 

Cnu 

2.411 

2.377 

2.343 

2.311 

2.277 

2.246 

Table  13.  Simulation  Results  for  Cno  and  Cdo  (AP09). 


Simulation  Data  (AP09) 

Mach 

0.5 

an 

0 

10 

20 

30 

40 

50 

60 

70 

80 

85 

Cl 

0.000 

0.268 

0.584 

1.462 

1.997 

2.449 

2.655 

2.464 

1.590 

0.828 

Cd 

0.097 

0.216 

0.406 

0.922 

1.421 

2.056 

2.793 

3.656 

4.377 

4.569 

Table  14.  Simulation  Results  for  Cl  and  Cd  at  High  AOA  (AP09). 
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